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Enhancement of critical current density in superconducting/magnetic multi-layers with slow 
magnetic relaxation dynamics and large magnetic susceptibility 
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We propose to use superconductor-magnet multi-layer structure to achieve high critical current density by 
invoking polaronic mechanism of pinning. The magnetic layers should have large magnetic susceptibility to 
enhance the coupling between vortices and magnetization in magnetic layers. The relaxation of the magnetiza- 
tion should be slow. When the velocity of vortices is low, they are dressed by nonuniform magnetization and 
move as polarons. In this case, the viscosity of vortices proportional to the magnetic relaxation time is enhanced 
significantly. As velocity increases, the polarons dissociate and the viscosity drops to the usual Bardeen-Stephen 
one, resulting in a jump in the I-V curve. Experimentally the jump shows up as a depinning transition and the 
corresponding current at the jump is the depinning current. For Nb and proper magnet multi-layer structure, we 
estimate the critical current density J c ~ 10 9 A/m 2 at magnetic field B as 1 T. 

PACS numbers: 74.25.Wx, 74.25.Sv, 74.25.Ha, 74.78.Fk 



I. INTRODUCTION 

One fascinating property of superconductors is ability to 
carry dissipationless current. With a transport current, vortices 
are induced inside the superconductor due to the magnetic 
field generated by the current. These vortices are driven by the 
Lorentz force exerted by the current and their motion causes 
voltage and dissipation. In inhomogeneous superconductors, 
the Lorentz force can be balanced by the pinning force due to 
defects. The strength of pinning thus determines how much 
dissipationless current the superconductor can carry, which is 
defined as the critical current. The random distributed point- 
like defects where superconductivity are weakened, are im- 
portant for pinning of vortices. 1 One may also introduce arti- 
ficially columnar pinning centers by heavy ion irradiation!^ In 
these cases, the pinning is caused by suppression of supercon- 
ductivity. 

An alternative approach to introduce pinning is to use mag- 
netic moments, which interact strongly with vortex. Such op- 
tion may be present in magnetic superconductorsP^The mag- 
netic moments can also be introduced artificially in hybrid 
systems consisting of superconducting and magnetic layers 5 . 
It was proposed in Ref. |6]that the domain walls can provide 
strong pinning with characteristic pinning energy OoM^. Here 
Oo = hc/(2e) is the quantum flux and is the magnetization 
at the wall. There are experimental attempts to enhance the 
critical current by putting magnetic particles^, dotd^ or fer- 
romagnet with domain walls on top of superconductors!^ 

Reduction of dissipation can be also achieved by enhance- 
ment of the vortex viscosity. At a given current J, the dis- 
sipation power for a superconductor without pinning due to 
quenched disorder is proportional to J 2 /?] where rj is the vor- 
tex viscosity. In nonmagnetic superconductors, rj is just the 
standard Bardeen-Stephen (BS) drag coefficient accounting 
for the dissipation in the normal vortex core. If one can in- 
crease significantly the vortex viscosity, superconductors can 
carry large current density with low dissipation, despite vor- 
tices are not pinned. It was shown that in magnetic supercon- 
ductors, motion of vortex lattice excites magnons^M When the 
kinematic condition £2(G) = G ■ v is satisfied, Cherenkov ra- 



diation of magnon occurs and the vortex viscosity is enhanced 
due to transferring energy into the magnetic subsystem, where 
energy is finally dissipated through magnetic damping. Here 
G is the lattice wave vector, v is the velocity of vortex lattice 
and £2(G) is the magnon spectrum. When the magnetic damp- 
ing is weak, magnetic domain walls are created dynamically 
due to the parametric instability and the viscosity is increased 
furtherpl 

Recently a polaronic mechanism of vortex pinning is pro- 
posed in Ref. [13]to explain the increase of critical current ob- 
served in ErNi2B2C below the incommensurate to commensu- 
rate spin density wave (SDW) transition at 2.3 KPl The tran- 
sition into the commensurate SDW phase leaves 1/20 spins 
free from molecular field 15 . These spins can be easily po- 
larized by vortices. These spins are Ising spins and experi- 
ence large crystal field splitting 16 , which results in slow re- 
laxation dynamics 17 . When the velocity of vortex lattice is 
low, a/v » t, the nonuniform component of free-spin magne- 
tization induced by vortex lattice follows the vortex motion, 
and the nonuniform magnetization and vortex form a polaron. 
Here a is the vortex lattice constant and r is the relaxation time 
for magnetizations. The effective viscosity of vortex lattice in- 
creases with the relaxation time. For a large velocity, a/v <K r, 
the nonuniform magnetization cannot follow the motion of 
vortex lattice and they are decoupled from each other. The vis- 
cosity of the system recovers to the conventional BS one. The 
decoupling or dissociation of polaron experimentally shows 
up as a depinning transition. The maximal critical current for 
ErNi 2 B 2 C is estimated as 10 10 A/m 2 at magnetic field B » 0.1 
T. The polaronic mechanism is also at work in other borocar- 
bides, cuprate and iron-based superconductors with magnetic 
rear earth ions locating between superconducting layers. 

The polaronic mechanism of pinning provides an additional 
routine to achieve high critical current. To optimize such pin- 
ning mechanism, we propose to use a multi-layer structure 
consisting of superconducting (S) and magnetic (M) layers 
shown in Fig. [T] to achieve high critical current. For that 
the magnetic layers should have high magnetic susceptibility 
at working magnetic field to ensure a strong coupling between 
magnetic moments and vortices. Secondly, the relaxation time 
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FIG. 1 . (color online) Schematic view of multi-layer structure con- 
sisting of alternating magnetic (M) layers (green) with thickness d,„ 
and superconducting (S) layers (blue) with thickness d s . The distri- 
bution of the magnetic field is shown by red lines. 



of the magnetization should be long. Thirdly, the penetration 
depth of the superconducting layers should be small. 



II. MODEL AND RESULTS 

Under external magnetic fields, the vortex lattice is induced 
inside the S layers. With a transport current, vortex lattice 
moves in response to the Lorentz force. In the quasistatic ap- 
proximation, the motion of vortex lattice is given by 



l 2 VxVxB+B = (Do^^r- r,-(f)] z, 



(1) 



where r,(f) = ro — vt is the vortex coordinate, z is the unit vec- 
tor along the z axis and A is the London penetration depth. In 
the flux flow region, the quenched disorder is average d out 
by vortex motion and the lattice ordering is improve d 18 * 19 ! 
The magnetic field inside the M layers is determined by the 
Maxwell equations 



V x (B - 4ttM) = 0, V • B = 0. 



(2) 



The magnetization M depends on B and is determined by the 
material properties. With a strong field and in static case, M 
is a nonlinear function of B and generally can be expressed as 
M(r) = Jiir 3 /(r-r',B(r')). The characteristic length of mag- 
netic subsystem is much smaller than A and we use a local ap- 
proximation /(r-r', B(r')) = <5(r-r')/(B(r')). B(r) has com- 
ponent uniform in space, Bo, and the other component nonuni- 
form in space, B(r), with Bo <K B. Thus the spatially nonuni- 
form magnetization M(r) is M(r) « 5/(B )/5fi B(r) = 
^o(Bo)B(r). In the following we assume the magnetic subsys- 
tem is isotropic and is characterized by a susceptibility xo (Bo) 
at Bo in static case. The magnetic field inside the M layer is 



determined by the equation V 2 B = 0. Since only the spatially 
nonuniform component M and B are responsible for pinning, 
we will focus on the nonuniform components in the following 
calculations. At the interface between the M and S layers, we 
use the standard boundary condition for the field parallel to 
the z-axis B z and field parallel to the interface B" 

B : \ S =B : \ M , B | l|s = (l-4^ ro )Bll| M . (3) 

Then we can obtain the magnetic field inside the M layers 

$0 exp (-iG x v x t) 



BUG>0,z) = a [e c z' + e- c ^]. [ + ^ . 
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with z' - z - n(d s + d,„), k s = V/L 2 + G 2 , and x' = (1 - 
4nx)~ 1 k s /G. Here n is the layer index and the vortex motion 
is assumed to be along the x direction. We consider square 
lattice G = (m x 2n/a,m y 2n/a) with a = V^oABo the lattice 
constant and m x , m y integers. 

We assume a relaxational dynamics for the M layers, 
M(cj) = ^(w)B m (w), with a dynamic susceptibility 



Xo 



1 + ia)T 



(6) 



Here we have assumed that the relaxation dynamics is gov- 
erned by a single relaxation time. This assumption is not es- 
sential but just for convenience of calculations. In the steady 
state, we have 



M(G 



,Z,t)= exp[(f - f)/-r] dt . 

Jo t 



(7) 



Because of the relaxation, M depends on the history of vortex 
motion. Due to slow relaxation of the magnetization, there is 
retardation between the time variation of induced nonuniform 
magnetization and vortex motion. As a result, the magneti- 
zation exerts a drag force to the vortex which is opposite to 
the driving force. The pinning force acting on a single vortex 
due to the induced magnetization in one M layer is given by 

F p - d ra J dxdy J d dzM ■ B m , which yields 



F p = 2^[l- exp (-2Gd m )] 



2a 2 xo 



GvtO> 2 



{\+A 2 G 2 fa 2 1 +(Gvt) 2 



(8) 



The I-V curve is determined by the equation of motion for 
vortex d s T]Bsv = d s Fi - F p with the electric field E = Bv/c 
and the Lorentz force Fi = 70o/c. Here 77^ is the BS vis- 
cosity tjbs = Q^IQ-n^^Pn) with p„ the resistivity just above 
T c and £ the coherence length. We consider a realistic case 
where aji7.ii) <K d m , d s . Taking into account only the domi- 
nant contribution G x — 2n/a and G y = in the summation, 
we obtain 



1+u 2 



(9) 
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FIG. 2. (color online) Calculated I-V curves for T p = 20 and T p = 2. 
For T p = 20 the system shows hysteresis in the I-V while for T p = 2 
no hysteresis is present. The green dotted line denotes the unstable 
solution. 
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with u = v/vq and vo = o/(2wt). 

At a small velocity u <K 1, the velocity is given by m = 
^l/(1 + ?>) which becomes inversely proportional to r for a 
large r. For a large u » 1, we recover the conventional BS 
viscosity v = Tl- The dependence of u on IFj, is shown in Fig. 
[2] Hysteresis is developped when J r p >%. For typical param- 
eters for Nb superconductor £ m A m 40 nm, p„ a 10~ 6 £2 ■ m 
and a — 40 nm at B » IT and xo = 0.05, T p > 8 requires 
t > 1 ps. For the relaxation time of order r » 1 jus, the ef- 
fective viscosity is enhanced by a factor of 10 6 compared to 
the bare BS one at v < a/r. Upon increasing the current, 
the velocity increases and at a critical current (velocity), the 
system jumps at dTildu — to the conventional BS branch 
due to the dissociation of vortex polaron. The effective critical 
current density for the whole system is given by 



1 



2- 4nx j (2nfAU s + d m 



(11) 



For d s — d m — 100 nm , we obtain J c a 10 9 A/m 2 . On the 
other hand, when one reduces the current from the conven- 
tional BS branch, the system jumps to the branches with high 
viscosity due to the formation of vortex polaron. We call this 
retrapping transition and the retrapping current J r is 



J r 



1 



1 - 27r;t"o 



Xamsads ac 1 
7tt A 1 An 2 d s + d m 



(12) 



For the parameters used before and t = 1 fj.s, we estimate 
J r « 2 x 10 6 A/m 2 . 



III. DISCUSSIONS 

Now we discuss the optimal thickness of M and S layers. 
For d m » a, the magnetic induction and the magnetization 
is practically uniform in the lateral direction in the middle of 
the M layer. This can be seen from Eqs. Q and Q that 
B m (G > 0) w exp(-2nd m / a) when -d m <K z' <s 0. As a 
result, the pinning force becomes practically d m independent 
when d m » a. In other words, the pinning is effective only 
near the boundaries between S and M layers in the area of 
thickness of the order a. On the other hand, the Lorentz force 
is proportional d s . Thus the effective critical current of the 
whole system J c is proportional to l/(d s + d m ) as described by 
Eq. 



(11 



Therefore the thinner of both M and S layers, the 
higher is the critical current of the system. 

Let us discuss the possible choice of S and M layers. The 
critical current decreases as A~ 4 because the smaller A, the 
more nonuniform is the magnetic field distribution inside the 
M layers, hence stronger pinning. Thus superconductors with 
smaller A are preferred. The critical current does not depend 
on r for sufficiently large r, while the viscosity in the branch 
with vortex polaron is proportional to t. The slow magnetic 
dynamics can be realized in spin glasses. Their relaxation 
is described by a broad spectrum of time scale, with aver- 
age time of the order 0.1 /zP^M For CuMn .o8, Xo ~ 0.002 
at B = 1 Tp2 One may enhance xq by tuning the concen- 
tration of magnetic metal in alloys. 23 One may use super- 
paramagnets with r as large as 1 s and with huge yp due to 
large magnetic moments in superparamagnetspS^l One may 
also use the recently synthesized cobalt-based and rare-earth- 
based single chain magnets with xo ~ 0.05 at B — IT and 
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Next we discuss the effect of quenched disorder. In the 
presence of quenched disorder, the vortex lines adjust them- 
selves to take the advantage of the pinning potential, which 
destroys the long-range lattice order. Below a threshold cur- 
rent, vortices remain pinned (actually they creep between pin- 
ning centers due to fluctuations). In this region, the polaronic 
mechanism does not play a role. When the current is high 
enough to depin the vortices from quenched disorder, vortices 
start to move and the lattice ordering is enhanced. By for- 
mation of polaron with the nonuniformly induced magnetiza- 
tion, the vortex viscosity is enhanced. At a critical velocity 
(current), the polaron dissociates and the system jumps to the 
conventional BS branch. Pinning due to quenched disorder 
works in the static region and polaronic pinning works in the 
dynamic region. The critical current of the whole system is 
the sum of these two threshold currents. Note that magne- 
tostriction in combination with quenched disorder enhance the 
polaronic pinning mechanism. 

The M/S multi-layer structure is naturally in some su- 
perconducting single crystals, such as RuS^GdC^O^ and 
(RE)Ba2Cu30 J 32 l 33 l where RE is the rear earth magnetic 
ions. In RuS^GdC^Og the magnetic moments order fer- 
romagnetically above T c thus the dominant enhancement of 
vortex viscosity is due to the radiation of magnons^. For 
(RE)Ba2Cu307, magnetic RE ions positioned between su- 
perconducting layers interact weakly with superconducting 
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electrons and order at very low Neel temperatures of the 
order 7# ~ IK. The polaronic mechanism is important 
above the magnetic ordering temperature, where spins are 
free. The London penetration depth of cuprate supercon- 
ductors is large A as 200 nm, thus the critical current is re- 
duced significantly compared to that for Nb multi-layer struc- 
ture, because J c drops as I /A 4 . Another natural realization is 
the recently discovered iron-based superconductors, such as 
(RE)FeAsOi_ x F x , where RE ions ordered antiferromagneti- 
cally below T N ~ 1KP 

IV. CONCLUSION 

To summarize, we have proposed superconductor-magnet 
multi-layer structure to achieve high critical current density 



based on the polaronic pinning mechanism. The critical cur- 
rent is estimated to be 10 9 A/m 2 at B » IT for an optimal 
configurations of Nb and proper magnet multi-layer structure. 
In the presence of quenched disorder, the polaronic pinning 
starts to work when vortices depin from quenched potential. 
Thus the total critical current of the system is the sum of de- 
pinning current due to quenched disorder and depinning cur- 
rent due to the polaronic mechanism. 
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